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Abstract 

We study the dynamics of systems with different time scales, when access only 
to the slow variables is allowed. We use the concept of Finite Size Lyapunov 
Exponent (FSLE) and consider both the case when the equations of motion for 
the slow components are known, and the situation when a scalar time series 
of one of the slow variables has been measured. A discussion on the effects 
of parameterizing the fast dynamics is given. We show that, although the 
computation of the largest Lyapunov exponent can be practically infeasible 
in complex dynamical systems, the computation of the FSLE allows to extract 
information on the characteristic time and on the predictability of the large- 
scale, slow-time dynamics even with moderate statistics and unresolved small 
scales. 
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I. INTRODUCTION 



In the last two decades, the problem of extracting information from a measured time 
series has been studied extensively, see e.g. ||lHT2|. Several attempts have been devoted to the 
issue of distinguishing between deterministic and stochastic behavior, where "deterministic" 
has to be interpreted as "dominated by a small number of excited modes" and "stochastic" 
as "dominated by a large number of excited degrees of freedom" . Once assessed the presence 
of low-dimensional chaotic dynamics, various methods have been devised for determining the 
statistical properties of the attractor and to build appropriate models for either predicting 
or describing the system evolution. 

Most methods for determining dynamical properties from measured signals are based on a 
procedure of phase-space reconstruction. Following the work of Packard et al and Takens 
0, the so-called time-embedding techniques have been developed to address this problem. 
Their use (e.g. via delay coordinates) allows, at least in principle, the determination of the 
dimensions 0, the Kolmogorov-Sinai entropy |^ and the Lyapunov exponents [0,11 of the 
system by the analysis of a time series of just one scalar variable. 

Unfortunately, this approach may have severe limitations in many practical situations. 
For example, the length of the time series is a crucial point in order to obtain reliable 
estimates of the phase-space properties of the system |]TT],|T2|. Further, there are simple 
stochastic processes that mimic a "false positive" answer to the search for low-dimensional 
chaotic dynamics, providing a finite value of the dimension under time-embedding in most 
practical cases [|T^ -|l5|. Analogously, simple systems characterized by on/off intermittency 
require additional care in the procedure of phase-space reconstruction and analysis 

Another problem is encountered in systems with many different time scales. In this case, 
it has been shown [ll8| , p!9| that the Lyapunov exponents may have a rather marginal role. 
The growth of a non-infinitesimal perturbation is indeed ruled by a non-linear mechanism 
which depends on the details of the system. For this reason, despite the positiveness of the 
largest Lyapunov exponent, it is possible to have a long predictability time for some specific 
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degrees of freedom. A typical example of this type of behavior is provided by 3D turbulence, 
that is characterized by the contemporary presence of a hierarchy of eddy turnover times. 
In this case, large-scale motions have a predictability time that is much larger than the 
one suggested by the value of the largest Lyapunov exponent. In such a situation, the 
predictability time for realistic perturbations may thus have no relationship with the growth 
rate of infinitesimal perturbations. 

As an attempt to overcome this problem, the concept of maximum Lyapunov exponent 
has recently been generalized in [l^ to the case of non-infinitesimal perturbations, intro- 



ducing the notion of the Finite Size Lyapunov Exponent (FSLE). In this work, we further 
elaborate on this issue and apply this method to the detection of "large scale" ("slow") 
dynamical properties of measured systems characterized by the contemporary presence of 
different time scales. In particular, we are concerned with systems that can be separated 
into a slow part, S, described by the phase-space variables Xs, and a fast part F, described 
by the variables Xf. The two subsystems are coupled through a term of typical strength e. 
In the limit e — 0, each of the two subsystems evolves independently with its own (chaotic) 
dynamics. The Lyapunov exponents of the slow and fast subsystem are Xg < Xf respectively. 

As for the coupling, we can either have a situation where the fast subsystem drives 
the slow one without being influenced by the latter, see for example [|l^,|l^, or a more 
generic coupling between the two parts [p^Sj^ . For the specific application we are concerned 
with here, the form of the coupling is not very important. Preliminary results on the 
predictability of a slow system S coupled with a faster system F have been discussed in the 
case of two coupled Lorenz models [^. In that work, the dynamics of the fast system was 



supposed to be known with arbitrary accuracy, and it was found that even if the value of 
the Lyapunov exponent is determined by the fast dynamics, the predictability of the slow 
system is dominated by its own characteristic time and it is almost unaffected by a small 
coupling with the fast dynamics. 

Physically, we may think of the fast subsystem as representing small scales that, both 
in real experiments and numerical simulations, are not resolved. Consistent with this inter- 



pretation, here we assume that the dynamics of the fast subsystem is poorly known, and 
investigate the effects of parameterizing the fast dynamics when one has access only to the 
slow dynamics. In this framework, we consider two different situations. In the first case, 
the equations of motion of the slow system are given. In the second case, we consider the 
computation of the FSLE directly from the a measured time series. The study of systems 
with two characteristic time scales is the first necessary step for the understanding of more 
realistic systems with several scales |pT |. 

The remainder of this paper is organized as follows. In Section 2 we discuss the notion 



of Finite Size Lyapunov Exponent introduced in |jT9[. In Section 3 we study the case of two 
coupled systems having different time scales, when access to the whole phase space of the 
slow system is allowed. In Section 4 we consider the same cases, but when just one scalar 
time series is supposed to have been measured. Section 5 gives conclusions and perspectives. 



II. EXTENSION OF THE LYAPUNOV EXPONENT TO FINITE 

PERTURBATIONS 

Here we recall the basic ingredients and the algorithm for computing the Finite Size 
Lyapunov Exponent (FSLE), referring to [jT9[ for further details. The definition of FSLE 
follows from that of error growing time Tr{6) for a perturbation of size 6. By definition, Tr{6) 
is the time that a perturbation with initial size 6 takes to grow by a factor r during the 
system evolution. In general, the perturbation with size 6 is supposed to be already aligned 
with the most unstable direction. The error ratio r should not be taken too large, in order 
to avoid the growth through different scales. In many applications, r = 2, so sometimes 
the Tr is also called the error doubling time. The Finite Size Lyapunov Exponent is defined 
from an ensemble average of predictability time according to 

where {■■■)t denotes the natural measure along the trajectory and (...) is the average over 
many realizations. The second equality comes from the definition of the time average along 
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a trajectory for a generic quantity A: 

{A), = U"Ait)dt = ^^ = ^. (2) 

in the particular case of A = 1/r [p!9|]. 

In the hmit of infinitesimal perturbations, 6 —>■ 0, this definition reduces to that of the 
leading Lyapunov exponent Amax- In practice, A (5) displays a plateau at the value Amax for 
sufficiently small S. 

In most systems, the smaller scales evolve faster, as in the classic example of three 
dimensional turbulent flows, and dominate the error growth for infinitesimal perturbations. 
When the size 6 of the perturbation cannot be considered any longer infinitesimal, all the 
scales whose typical size is smaller than 6 experience a diffusive separation and do not 
contribute to the exponential divergence in phase space. At this stage, the behavior of X{6) 
is governed by the nonlinear evolution of the perturbation, and, in general, A (5) < Amax- The 
decrease of A (5) does follow a system-dependent law. In some cases, A (5) can be predicted 
by dimensional considerations. For the fully developed turbulence, for example, dimensional 
considerations lead to the universal law X{6) ~ in the inertial range [|1^ . 

Therefore, the behavior of A as a function of 6 contains important informations on the 
characteristic times governing the system, and it is a powerful tool for investigating the 
behavior of high-dimensional dynamical systems involving many characteristic scales in space 
and time. 

To practically compute the FSLE, one has first to define a series of thresholds 5„ = r"5o, 
and to measure the time Tr((5„) that a perturbation with size 6n takes to grow up to Sn+i- 
The time Tr(5„) is obtained by following the evolution of the perturbation from its initial 
size 5min up to the largest threshold 5max- This can be done, for example, by integrating two 
trajectories of the system that start at an initial distance S^in- In general, one must take 
<^min <^ ^0; ordcr to allow the direction of the initial perturbation to align with the most 
unstable direction in the phase-space. The FSLE, A((5„), is then computed by averaging the 
predictability times over several realizations, see equation (|l|). 



Note that the FSLE has conceptual similarities with the e-entropy |2^. This latter 
measures the bandwidth that is necessary for reproducing the trajectory of a system within 
a finite accuracy 6. The e-entropy approach has already been applied to the analysis of 
simple systems and experimental data |^ , giving interesting results. The direct calculation 
of the e-entropy, however, is much more expensive than that of the FSLE. This latter, in 
fact, is not more expensive than that of the largest Lyapunov exponent Amax- 



III. FSLE AND SMALL SCALE PARAMETERIZATION 

In this section we study the case of a slow system S, described by the variables Xs, coupled 
with a fast system F described by the variables x j . The equations of motion governing the 
slow variables are supposed to be known, and we study the effects of parameterizing the fast 
dynamics. 

To study the evolution of the perturbation, we consider two trajectories x = {Xs,Xf) 
(reference) and x' (perturbed) starting from two nearby locations in phase space. The 
perturbed trajectory is then made to evolve according either to the same equations as the 
reference one, or to modified equations where the fast dynamics is replaced by a stochastic 
process, or simply neglected. 



A. Coupled maps 

The first example is provided by two coupled maps, namely 

Xs{n + 1) = {1 - e) fs[xs{n)] + eg[Xs{n),Xf{n)], modi 
Xf{n + 1) = {1 - e) ff[xf{n)] + eg[xf{n),Xs{n)], modi 
where fs and // are maps of the unit interval [0, 1] into itself. Here we use 



(3) 



fs i^s 



e^' Xs 



mod 1 



ff{xf) = e'f Xf 



mod 1 



g{xs,Xf) = cos{27i{xs + Xf)) (4) 



with As < A/. Equations (^) completely define the dynamics of the system. We assume, 
however, that we can have access only to the slow variable Xs{n). The FSLE has thus to be 
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computed only from the time evolution of Xs{n). 

Before discussing how this can be achieved, we note that equation is inadequate in 
the case of maps. Definition (|1|), in fact, tacitly assumes that we are able to determine the 
time when the size of the perturbation is exactly equal to the fixed threshold 6. In the case 
of maps, this may not be possible. The appropriate definition in this case thus becomes (see 

H) 

where 6 is the initial size of the perturbation and 6{nr) is its size at the (discrete) time rir- 
Here rir is the time at which the size of the perturbation first gets larger than (or equal 
to) r6, i.e., 6{nr) > r6 and 6(71^ — 1) < r6. The average (■ ■ ■) is over an ensemble of many 
realizations, as in the original definition (|l|). 

Let us now discuss how to compute X{6) from the knowledge of Xs{n). From a point 
{xs,Xf) on the system's attractor, we generate a new point representing the perturbed 
trajectory (5min < 1) 

Xg = Xg -\- 5min; X j: = Xf + ^rnin (6) 

and iterate the coupled maps for the original trajectory and the perturbed one. Note that, 
in this case, the perturbation has been applied to both the slow and the fast variables. We 
then compute A (5) from eq. (^. 

Figure 1 shows the value of X{6) versus 6 for the coupled maps system. The curve 
with filled triangles has been obtained by defining the distance in phase space as 5 = 
\xg{n) — x'g{n)\. The curve denoted by the filled squares has been obtained with a different 
definition of distance, namely 5 = {[xg{n)—x'g{n)]'^ + [xg{n — l)—x'g{n — l)]'^}^^'^ (reminiscent 
of the time embedding procedure, see next section). Both curves are obtained by an average 
over 10^ samples for each value of 6; analogous results are obtained with more limited 
statistics. For small S, the dynamics of the perturbation is driven by the fast mode, and 
X{6) tends toward Amax — ^/ = 0.5. For large values of 6, the growth of the perturbation 



is governed mainly by the slow dynamics and X{S) approaches = 0.1. The transition 
between the two regimes takes place at 5 ~ e = 2 x 10^'^; changing the strength of the 
coupling modifies only the value of 6 where the transition takes place. It is worth noting 
that the two definitions of distance used in obtaining figure 1 give almost coincident results. 
This indicates that, at least in this case, the definition of distance which is employed is not 
crucial, at variance with what happens in other cases (e.g. for on/off intermittent systems 
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As discussed above, we have always used (5min ^ in order to allow the direction of the 
initial perturbation to align with the most unstable phase-space direction, and hence to be 
able to recover the largest Lyapunov exponent in the limit of infinitesimal perturbations. In 
figure 2 (curve with filled triangles) we show what happens when this prescription is relaxed. 
In this case, the value of A (5) is underestimated for 6 ~ 6ram- This is due to the fact that 
the perturbation is not along the most unstable direction. The process of alignment of the 
perturbation along the unstable direction may result in a decrease of the distance between 
the two trajectories at initial times (i.e., at small 6), leading to a lower value of A(5). The 
computation of X{6) without relaxation can indeed be more appropriate for characterizing 
short time predictability with large initial error, but we cannot expect, in this case, to 
asymptotically recover the largest Lyapunov exponent. 

From these results, it is apparent that the effects of the fastest dynamics are seen only 
when the size of the perturbation is small enough. As a consequence, one has that the pre- 
dictability for finite-size perturbations may be unaffected by the particular parameterization 
of the faster scales. To confirm this inference, in figure 3 we show X{6) for a case where 
in the evolution of the perturbed trajectory (x(,,Xj) the fast variable is replaced by a 
sequence of random numbers uniformly distributed in the interval [0, 1] (curve with filled 
triangles). The fact that x'j is now a random variable can be detected only for small enough 
6 (where X{6) is ruled by a logarithmic law, see []T^|). The characterization of the large scale 



dynamics, however, is unaffected by the incorrect parameterization of the fast dynamics. 
Note, also, an interesting point. The results in figure 1 indicate that the small-scale 
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dynamics is correctly recovered when a one-dimensional time series is used (xs(n)), even 
though the full dynamics is two-dimensional {xs{n),Xf{n)). This event, possibly surprising 
at first sight, is due to the fact that the perturbed trajectory has been obtained by acting 
on the full phase space of the system; i.e., the perturbation has been made on both Xs{n) 
and Xf{n). When only Xs{n) is perturbed, as in the curve shown in figure 3, it is not 
possible to recover the fast small-scale dynamics without resorting to methods such as the 
time embedding technique. In this case, however, other problems appear, as discussed in 
the next section. 



B. Coupled Lorenz models 



To illustrate the application of the FSLE technique to the case of continuous-time dy- 



namical systems, here we consider a system obtained by coupling two Lorenz models 
having time scales that differ by a factor a. 

The slow subsystem is coupled through the Rayleigh number R to the fast one; for 
simplicity, the fast subsystem does not feel any feedback from the slow one. More generic 



(small) couplings do not qualitatively change the results, see |2^. The equations for the 
whole system are 



dt 



-axs + ays 



(7) 



^ = -XsVs + (-R + (^Zf)Xs - Us 

-^sys bzg 

^ = {-(TXf + ayf)-a 

= i-^fyf + R^f-yf)-(^ 

^ = ixfyf-bzf)-a 

where the parameter e controls the strength of the coupling and a the relative time scale. 

In the first type of simulations, both the reference and the perturbed trajectories evolve 
with the same equations of motion (|^. Again, we assume that only the slow variables are 
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accessible; the norm is defined as the Euchdean distance in the three dimensional space 

The parameters used in the numerical integrations are R = 45.92, a = 16, b = A; the 
Lyapunov exponent of the slow subsystems is ^ 1.5. The ratio of the time scales of 
the two systems is a = 5, hence Xj — Amax — 7.5. The results for e = 10^"^ are shown 
in figure 4. As in the case of the coupled maps, A(5) displays two plateaus at \{6) — A/ 
and X{S) ~ As, corresponding respectively to the fast and slow dynamics, and a transition 
region at 5 ~ e. Again, the small-scale dynamics (associated with the fast variables) can 
be recovered because we have perturbed the trajectory in the full, six-dimensional, phase 
space. 

To investigate the role of small scale parameterization, also for this system we have 
considered a situation where the "true" dynamics of the fast variables x'j,y'j,z'j in the 
perturbation is replaced by a stochastic process, i.e. 

dx' I I / 

^ = -x',y', + {R + ev)x', - y', (8) 



dt 

dz^ 
dt 

where rj is a Gaussian white noise process with variance equal to that of the fast component 
in the original system (|^). Analogously, we have considered a case where the fast dynamics 
of the perturbation is simply neglected; this corresponds to taking e = in the evolution 
equations for the perturbation x'^, y'^, z'^. The two corresponding curves of X{6) are shown in 
figure 4 (filled triangles and open diamonds). 

The results shown in figure 4 confirm that the estimate of X{6) for the slow variables is 
practically unaffected, for large values of 6, by the details of the fast dynamics. A similar 
result was obtained in for the chaotic or stochastic resonance of a driven nonlinear os- 



cillator. In the present situation, one may even neglect the fast dynamics, and still obtain a 
reliable estimate of the slow evolution (and of the Lyapunov exponent associated with the 
slow variables). In particular, the lack of knowledge of the fast dynamics has an effect which 
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is similar to that associated with the presence of noise. The inadequacy of the parameter- 
ization for scales smaller than e is is reflected, in both cases, in a large value of the FSLE, 
i.e., in a poor predictability of the phase-space dynamics on small scales. At larger scales 
and slower times, the FSLE coincides with the Lyapunov exponent of the (uncoupled) slow 
subsystem. This defines the predictive skill of the "incomplete model" on those scales. 



IV. COMPUTATION OF THE FSLE FROM MEASURED DATA 

In the case of measured time series, it is not usually possible to have access to the whole 
set of variables describing the system. Consistent with these limitations, here we suppose 
that only one time series of a scalar observable quantity h^, function of the slow phase-space 
variables of the system, is given. Additionally, in most experimental situation, the time 
series of /i„ is characterized by limited statistics. 

The first step is the procedure of phase-space reconstruction. The time-embedding 
method allows to reconstruct a pseudo phase space with dimension M, by using time 
delay coordinates of the observed variable. A vector in this phase space is then defined as: 

Xn = (hn, hn-T, • • • , ^n-(Af-l)T) (9) 

where r is a suitably chosen time delay, see e.g. p|J^ for a discussion on the optimal choice 
of T. The method for computing the experimental FSLE is then a simple modification of the 
standard algorithm for the Lyapunov exponent [T^ which measures the average separation 
between trajectories in the embedding space. 

For each reconstructed vector X„, its nearest neighbor Xm is determined. If the sepa- 
ration 6 = \Xn — Xm\ is smaller that a given threshold 5miru the trajectories starting from 
Xn and Xm are used to compute the FSLE, according to the algorithm discussed above. 
As in the case of maps, the trajectory is not continuous in time and one has to adopt the 
definition (^. 

Also in this case we require (5min to be considerable smaller than 60, to allow the vector 
separating the two trajectories to align with the maximally expanding direction. Clearly, 
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this may severely limit the available statistics. A trivial geometrical argument shows that 
the probability of finding two points at a given distance in the embedding space becomes 
extremely small for high embedding dimensions. We have found this lack of statistics to be 
the most important limitation that prevents from taking the limit 6^0, and thus from 
estimating the largest Lyapunov exponent. On the other hand, the statistics grows with the 
perturbation threshold 6j, we may thus expect to be able to compute \{Sj) for sufficiently 
large values of 6j. Another crucial point is related to the fact that the number of degrees 
of freedom which participate in the slow dynamics is (usually much) smaller than the total 
number of excited modes. Thus, the embedding dimension which is required for estimating 
X{6) for moderate values of 6 is smaller that than needed for estimating the largest Lyapunov 
exponent. 

As a first example, we consider the case where the signal /i„ is generated by the component 
Xs{n) of the coupled maps model (^. The coupling parameter is herein chosen to be e = 0.02 
(larger than that used in the previous section), in order to be able to study the small-scale 
behavior at 6 < e even with moderate statistics. The results are shown in figure 5. The 
three curves show the FSLE as obtained in the previous section (filled squares), and that 
from the time-embedding method with embedding dimensions M = 1 (filled triangles) and 
M = 2 (open diamonds) with time delay r = 1. As expected the computation of 

the Lyapunov exponent requires in this case an embedding dimension M = 2 in order 
to resolve the fast dynamics. In fact, at variance with the results shown in the previous 
section, the one- dimensional time series is not enough here because we have no control on 
the perturbation in the fast variable. It is interesting, however, that one can obtain X{6) for 
large S already with an embedding dimension which refiects the dimensionality of the slow 
system (M = 1). 

Note that, in order to have good statistics, each point if figure ^ has been obtained by 
averaging (|^) over 10^ samples. This requires a time series of about 10^ points for M = 1 
and more than 10^° points in the case M = 2 to resolve the small scales at 5 ~ 10~^. 
Thus, although it is in principle possible to extract information on the fast dynamics and 
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on the largest Lyapunov exponent from a measured time series, the statistics required is 
so prohibitive that this procedure may be infeasible in reahstic situations. On the other 
hand, it is possible to extract information on the large-scale Lyapunov exponent with an 
embedding dimension of the order of the number of degree of freedom involved in the slow 
dynamics. 

We observe that this result could not be obtained by simply neglecting the fast component 
as in a measured time series one has no direct access to the equations of motion. 

As a second example, we now apply the same machinery to the case of the two coupled 
Lorenz systems described by equation (^. Figure |^ shows the results obtained by using the 
variable Xg of (J^), with M = 3, r = 0.02 and a total number of = 500, 000 points in the 
time series. The coupling constant between the models is now e = 0.05. The perturbation 
threshold is fixed as 6min = 0.005 and 6o = 0.05. The plateau corresponding to the large 
scales is clearly visible, while, in spite of the large number of points, the contribution of the 
fast system is not resolved. Clearly, this would require a larger embedding dimension and 
an increase of the smallest resolved value of 6, at the cost of an unrealistic increase of the 
number of points in the time series. 

V. CONCLUSIONS 

In this work we have discussed how to use the method of the Finite Size Lyapunov Expo- 
nent (FSLE) for determining the slow dynamics of systems with many different characteristic 
times. In particular, we have considered the case when full access to the slow dynamics is 
allowed, and the more realistic case when just one scalar time series of a slow variable has 
been measured. The basic idea is to compute the FSLE, A (5), as a function of 6 in the 
framework of the embedding technique. In this case, the behavior of X{6) at large values of 
6 gives information on the Lyapunov exponents associated with the slow dynamics. 

By contrast, the behavior of X{6) for small value of 6 gives information on the fast 
dynamics. By considering the limit of A (5) for 5 ^ 0, it is possible, at least in principle, 
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to discriminate between "stochastic" systems and "chaotic but deterministic" ones. For an 
accurate estimate of X{6) at small 5's, however, it is necessary to use a very large number 
of points in the time series. This fact makes practically infeasible the calculation of the 
largest Lyapunov exponent (associated with the fast dynamics) in most complex dynamical 
systems. 

Nevertheless, very often the slow, large-scale dynamics is the most interesting one phys- 
ically. The results obtained here indicate that the slow dynamics may be satisfactorily 
detected even with a limited number of points and a moderate embedding dimension. In 
these systems, one can thus obtain a satisfactory prediction for the slow, physically inter- 
esting scales even when access to the (much more unpredictable) fast scales is not available. 
This also indicates that, at least in the examples considered here the parameterization of the 
fast time scales seems not to be crucial, as the internal dynamics of the slow modes plays 
the dominant role. 

One could wonder, then, how general the results presented in this paper are. Previous 
works on more complex theoretical models indicate that the crucial point is not the 
dimensionality of the system or the details of the couplings, but rather the existence of well 
separated, weakly interacting, scales. In the present work we have shown that the FSLE 
technique may be successfully applied also in the case of measured time series. 
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FIGURES 

FIG. 1. The FSLE, X{6), as a function of 6 for the coupled maps ^) with = 0.1, A/ = 0.5 
and e = 2 X 10"^. The two curves refer to the different definitions of distance discussed in the text. 
The parameters of the perturbation are (^min = 10~^, 6o = 10~^, 6max = 0.1, r = 2 and the average 
is over 10^ realizations for each point in the FSLE curve. The horizontal lines indicate the values 
of Xs and Xf. 

FIG. 2. The FSLE, X{6), as a function of 6 for the coupled maps ^) with the same 
parameters as in figure ||. The filled squares refer to the results obtained with (5min = 10~^ 
and So = 10~^ (same curve as in figure ||). The filled triangles show the behavior of X{6) when 
Smin = <5o = 10^^ and the perturbation did not start along the most expanding direction. 

FIG. 3. The FSLE, X{S), as a function of S for the coupled maps (^, |^. The parameters are 
as in figure 0. The filled squares are as in figure 1. The filled triangles refer to the case where the 
fast variables are replaced by a sequence of random numbers uniformly distributed in [0, 1]. 

FIG. 4. The FSLE, X{6), as a function of 6 for the coupled Lorenz systems with R = 45.92, 
a = 16, b = A and coupling e = 10^"^ (filled squares). The filled triangles indicate the results ob- 
tained when the perturbed trajectory evolves according to the modified dynamics (^) with e = 10~^. 
The open diamonds represent the case e = in the perturbed trajectory (no fast dynamics). 

FIG. 5. The FSLE, X{6), as a function of 6 for the time series x(n) obtained from the coupled 
maps (^, §) with ^min = 10~^, Sq = and coupling e = 0.02 (filled squares). The triangles 

indicate the results for an embedding dimension M = 1 and the diamonds for M = 2. The number 
of realizations used for each point in the FSLE is 10^. The total number of point in the time series 
is 10^ for the case M = 1 and more than 10^*^ for M = 2. 

FIG. 6. The FSLE, X{6), as a function of 6 for the time series Xs{t) obtained from the coupled 
Lorenz models (0) with coupling e = 0.05 (filled squares). The filled triangles indicate the results for 
an embedding dimension M = 3, time delay r = 0.02 and a total number of points = 500, 000. 
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